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@ Method for evaluating parameters related to the elastic properties of subsurface earth formations. 

(g) Three-axis vertical seismic profile measurements are 
analyzed to derive values for the velocity v p and angle of 
incidence of a plane compressional wave at the borehole 
and for the shear wave velocity v s and angle of incidence 4^ of a 
plane shear wave at the borehole. In the described embodiment 
this is accomplished by finding values for these parameters 
which maximize the expression 

E«Ka?*3*) H aM a 1lap*fiiF - |Ca p -a s )KuCco)Ma ? -a 3 ||2 

where 

0) is angular frequency; 

a p and a s are column vectors representing compressio- 
nal and shear wave particle motions at frequency o at a 
multiplicity of locations in the borehole; 

u(co) represents components in a plane of waveform 
measurements obtained at these locations; 

II i| represents the Euclidean norm of a complex vector; 
and 

() H denotes the complex conjugate transpose of a 
matrix. 
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D scription 



METHOD FOR EVALUATING PARAMETERS RELATED TO THE ELASTIC PROPERTIES OF SUBSURFACE EARTH 

FORMATIONS 



BACKGROUND OF THE INVENTION 



This invention relates to methods for evaluating 
parameters related to the elastic properties of 
subsurface earth formations, in particular the veloc- 
ity and direction of propagation of compressional (P) 
and shear (S) wave components of acoustic energy 
propagating through the formations. The invention 
also provides a way of separating the individual 
waveforms of such components. 

Seismic investigation of earth formations is a 
well-established means of deriving information 
about sub-surface geological conditions, and is a 
valuable tool in the search for natural resources 
such as oil, gas and ores. In such an investigation a 
seismic wave is propagated through the earth 
formation to be investigated, by means of a source 
such as an explosion or a mechanical vibrator. One 
or more acoustic detectors spaced from the source 
are used to sense the propagating seismic wave as it 
passes different locations and comparison of waves 
detected at these different locations enables infor- 
mation about the formation to be inferred. In one 
form of seismic investigation the source and/or the 
detectors are placed in deep boreholes drilled into 
the earth formation, enabling a vertical seismic 
profile (VSP) or well-to-wel! profile to be obtained. 
The detectors are typically 3-axis geophones which 
are placed in firm contact with the wall of the 
borehole and which can produce electrical signals 
indicative of movement of the formation along each 
of three orthogonal axes. 

The waveform signals thus derived can be ana- 
lyzed in a variety of ways to obtain information about 
subsurface conditions. Recently interest has been 
increasing in the possibility of deriving individual 
values for the compressional and shear propagation 
characteristics of formations. Whereas fluids can 
only support compressional (P) sonic waves, in 
which the particle motion comprising the wave is 
parallel to the direction of wave propagation, solids 
can also sustain shear (S) waves, in which the 
direction of particle motion is transverse to the 
direction of wave propagation. The velocities of P 
and S waves are two independent parameters 
related to the elastic properties of formations and 
are of interest for various reasons. Thus one useful 
indicator of lithology type is the ratio of P and S wave 
velocities. Imaging of reflected compressional and 
shear waves in surface seismics or vertical seismic 
profiles can provide valuable information, but re- 
quires values of P and S wave velocity as a function 
of depth. 

Techniques are already known for separating and 
extracting P wave components and S wave compo- 
nents in the waveforms obtained during VSP 
measurements. See N for example US Patent 



No. 4 648 039 to Devaney et al, assigned to the 
5 assignee of this application. However, these tech- 
niques require knowledge of the velocity of propaga- 
tion of the P and S waves in the formation. 
Determination of P wave velocity is possible by 
having a skilled user pick waveform events which 

10 appear characteristic of the arrival of P waves at the 
detectors. The time offset between corresponding 
arrivals at detectors having a known spacing is 
related to the velocity. 

However, picking arrivals accurately is difficult and 

15 time-consuming even in the case of P wave events, 
which are generally the first significant signals to 
appear on a seismic recording. S wave events 
usually occur later, in the midst of other events such 
as P wave reflections and mode conversions. Thus 

20 picking S wave arrivals can be very difficult and is 
sometimes impossible. 

Moreover, if the seismic energy is not propagating 
vertically, as is the case in VSP measurements with 
an offset source, it is necessary to determine the 

25 angle of incidence of the waves at the detectors as 
well as the time offset between between detectors in 
order to obtain correct values for the velocity of 
propagation. Although the angle of incidence can be 
determined to some extent from the particle motion 

30 of the waves within time windows specified by an 
experienced user, this operation again is difficult, 
time-consuming and subject to error, as well as 
requiring the intervention of a skilled user. 

Accordingly there is a need for a method of 

35 determining velocities of both P and S wave 
propagation which is fast and reliable, which re- 
duces or eliminates the need for user interaction and 
skill, and which avoids the need to identify specific 
events such as arrivals. 

40 It is an object of this invention to provide a method 

of determining velocity of P wave and S wave 
propagation in a formation from three (or two) 
component seismic measurements such as VSP 
measurements, which reduces or avoids the need 

45 for identification of arrivals. 

It is another object of this invention to provide a 
method of determining direction of P wave and S 
wave propagation in a formation from three (or two) 
component seismic measurements such as VSP 

50 measurements. 

It is a further object of this invention to provide a 
method of separating P and S waveform compo- 
nents in seismic measurements such as VSP 
measurements. 

55 

SUMMARY OF THE INVENTION 



60 According to one aspect of this invention there is 

provided a method of evaluating parameters related 
to the elastic properties of a subsurface earth 
formation traversed by a borehole. Measurements 
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are obtained along at least two, and preferably three, 
orthogonal axes from a seismic profile operation in 
the borehole. This profile operation may be, for 
example, a vertical seismic profile or an offset 
vertical seismic profile, and may be performed with 
the source on the surface or in the borehole. Seismic 
propagation through the formation is described in 
terms of a model incorporating a plane compressio- 
nal wave and a plane shear wave propagating with 
unspecified velocity and direction. Values are then 
derived for the velocity and direction of propagation 
of these plane waves that provide a best fit of the 
model to the measurements. Finally a tangible 
record of the derived values is produced. 

Typically the values for the velocity and direction 
of propagation are derived to provide a best fit of the 
model to the measurements obtained at a small 
number (for example, between three and eleven) of 
adjacent measurement locations. For typical spac- 
ings between measurement locations even a devi- 
ated borehole may be incorporated in the analysis, 
since the curvature of a typical deviated borehole 
over the distance spanned by such a number of 
locations is so slight that it can be ignored. 

It is believed that the inventor hereof is the first to 
suggest the possibility of analyzing seismic propa- 
gation in terms of the behavior of plane wave 
components over short sections of a borehole. This 
analysis is based on the assumptions that the 
measurements can be appropriately described in 
terms of propagation of plane wavefronts of seismic 
energy; that the formation is locally homogeneous 
and isotropic; and that there is no coupling between 
orthogonally polarized shear wave components. The 
inventor hereof has found that these assumptions 
are sufficiently well satisfied in most situations for 
the results produced by the method of the invention 
to be of practical benefit. 

The values of compressional wave velocity v p and 
angle of incidence *F P at the borehole and of shear 
wave velocity v s and angle of incidence *F$ at the 
borehole may be derived in accordance with the 
expression: 

maXvp/Kp v s .y s I«|(a p + a s ) H u(o))| 2 /||a p + a s || 2 + 
|(a p -a s ) H u(co)| 2 /||a p -a s || 2 
where co is angular frequency; 

ap and a s are column vectors representing 
compressional and shear wave particle motions at 
frequency o at a multiplicity of locations in the 
borehole; 

u(co) represents components in a plane of wave- 
form measurements obtained at the locations: 

II II represents the Euclidean norm of a complex 
vector; and 

() H denotes the complex conjugate transpose of a 
matrix. 

Thereafter^ separate compressional and shear 
waveforms u p (co) and u s (o>) may be derived in 
accordance with the expressions 

^up(co) = [Up + £s)/ll£ P + £ s || 2 4- (a p 
-i^/|U p -j_ 9 ||2]Hu( C o| ^ ~ 

^upfco) - l(7 p + a s )/ll a p + £ 0 || 2 -Up -a s )/||a p 
-iufi 2 ] H u(cD) ~ ~ " ~~ " ~ 

where ^ p and ^ s are values of the column vectors a p 
and a s corresponding to the values derived for the 



velocity and direction of propagation of said plane 
waves. 



BRIEF DESCRIPTION OF THE DRAWINGS 



Additional objects and features of the invention 
will become more apparent upon consideration of 
10 the following detailed description of the invention, 
reference being had to the accompanying drawings 
in which: 

Figure 1 is a schematic diagram of a borehole 
seismic operation for collecting data for use in 
15 accordance with this invention : 

Figure 2 is a schematic diagram showing a 
coordinate system used in an embodiment of 
the present invention: 

Figure 3 is a diagram showing directions of 
20 propagation and particle motion of P and Sv 

waves in the coordinate system of Figure 2; and 
Figures 4a and 4b show a flow diagram of a 
method in accordance with this invention for 
evaluating parameters related to elastic proper- 
25 ties of a formation . 



30 



DETAILED DESCRIPTION 



Referring to Figure 1, a logging tool 10 is 
suspended on an armored cable 12 in a borehole 14 
penetrating an earth formation 16. The borehole 14 is 
filled with liquid 18, for example drilling mud used to 

35 stabilize the borehole 14 (in the case of an open, 
uncased hole) and prevent formation fluids from 
escaping up the borehole. The cable 12 extends up 
the borehole, over a sheave wheel 20 in a derrick 21 
to a winch forming part of surface equipment 22. 

40 Known depth gauge apparatus (not shown) is 
provided to measure cable displacement over the 
sheave wheel 20 and thus depth of the logging tool 
10 in the borehole 14. 
The tool 10 may be designed and constructed in 

45 accordance with the teaching of US Patent 
4 563 757, which is incorporated herein by reference. 
In particular the tool has a retractable anchoring pad 
24 for urging the body of the tool 10 firmly against 
the wall of the borehole 14. The body of the tool 

50 includes a three-axis geophone system arranged to 
detect particle motion along each of three orthogo- 
nal axes and to produce electrical signals indicative 
thereof. The geophones are gimballed. as described 
in US Patent Application No. 855 024. so that two of 

55 them have horizontal axes and the third a vertical 
axis irrespective of orientation of the tool 10, and a 
device is included in the tool 10 to produce a signal 
indicative of orientation of the body of the tool 10. 
Processing and interface circuitry within the tool 10 

60 amplifies, samples and digitizes these signals for 
transmission and communicates them to the surface 
equipment 22 via the cable 12. Electrical power and 
control signals for coordinating operation of the tool 
10 are generated by the surface equipment 22 and 

65 likewise communicated via the cable 12 to appropri- 
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ate circuitry provided within the tool 10. 

A source of acoustic energy 26 is located on the 
surface, near the borehole 14. This source may be, 
for example, an explosive device or an airgun in a 
shallow liquid-filled hole, or a vibrating device 
pressed against the earth's surface, and is con- 
trolled by signals transmitted along a line 28 from the 
surface equipment 22. Details of suitable sources 
are well known to those skilled in the art and need 
not be described here. 

During acquisition of data, the logging tool 10 is 
positioned within the borehole 14, the pad 24 is 
extended to clamp the tool in place and the source 
26 is energized to propagate seismic energy into the 
formation 16. This energy reaches the borehole 14, 
both directly and after reflection from interfaces 
between different layers constituting the formation 
16. The seismic energy causes displacement of the 
particles comprising the formation 16, and this 
displacement is sensed by the geophones in the 
body of the tool 10 to generate electrical signals 
representative of the motion along each of the three 
orthogonal axes. These signals are transmitted via 
the cable 12 to the surface equipment 22 which 
receives and records them on chart and/or magnetic 
tape recorders as a function of time. In addition the 
equipment 22 may, as described below, analyze the 
data represented by these signals to yield further 
information which is also recorded. 

When recording of the geophone signals has 
terminated the pad 24 is released and the cable 12 is 
moved by the winch to bring the tool 10 to another 
depth in the borehole 14, and the procedure is 
repeated. The distance through which the cable 12 is 
displaced depends on the subsequent processing of 
the recorded signals, which typically requires signals 
to be recorded at points separated by no more than 
a certain maximum spacing. This spacing is depend- 
ent in known manner on the minimum wavelength of 
the seismic signals of interest. 

Further details of the operation of the tool 10 may 
be found in US Patent 4 563 757. Those skilled in the 
art will recognize that other ways of obtaining the 
required data may be used. Thus, for example, an 
array of geophones may be used instead of a single 
assembly. 

The surface equipment 22 typically incorporates a 
data processor 30 for coordinating and controlling 
the logging operation, and this processor may also 
be used for analysis of the recorded seismic signals 
at the wellsite. Alternatively or in addition, the 
recordings may be transferred to a remote location 
for subsequent more detailed analysis. It will be 
understood by those skilled in the art that this 
analysis can be implemented, for example, by 
appropriate programming of a general purpose 
digital computer or by means of special purpose 
electronic circuitry. 

The purpose of the analysis is primarily to estimate 
values for four parameters related to the elastic 
properties of the formation. These parameters are 
v p the velocity of compressional (P) waves in the 
formation; 

*F P the direction of propagation of P waves; 
v s the velocity of shear (S) waves in the formation ; 



H's the direction of propagation of S waves. 
The analysis is performed in terms of plane wave 
components propagating in a (generally vertical) 
plane containing the source 26 and the measure- 
5 ment locations of the too! 10 in the borehole 14. Thus 
as shown in Figure 2 the borehole 14 and the 
measurement locations 32 are considered to extend 
along an axis z, while the line joining the top of the 
borehole 14 to the position of the source 26 extends 

10 parallel to an axis x. The wave components whose 
parameters are determined are those components 
which propagate along directions contained in the 
x-z plane. Accordingly, the parameters v s and 4^ 
relate specifically to a shear wave component with 

15 particle motion [perpendicular to the propagation 
direction, that is to an Sv(vertically polarized) shear 
wave (see Fjgure 3). Values of the four parameters 
for the depth level of each measurement location are 
determined using the measurements for the (2M + 1 ) 

20 locations straddling that depth level. Deviated 
(non-rectilinear) boreholes can also be accommo- 
dated in this analysis since the curvature of a typical 
deviated borehole over the distance spanned by the 
(2M + 1) measurement locations is so slight that the 

25 interval can be considered as approximately straight. 
The parameter values are estimated in effect by 
finding values which produce a minimum error value 
E in the expression 

E = £J|u(co) -£(co)||2 (1) 

30 where 

u(co) represents the components in the x-z plane 
of the actual waveform measurements obtained with 
the tool 10 at the (2M-M) measurements locations 
32,^0) being angular frequency; 

35 iL(*°) represents model waveform components in 
the x-z plane derived from estimated values of v p , 
Tp, v s and Ts; and 

II II represents the Euclidean norm of a complex 
vector. 

40 The waveform components Gjco) are formally 
defined (although not actually derived) in terms of 
the frequency domain components u p (ca) and u s (o>) 
of separate P and Sv waves respectively at the 
center measurement location, as follows: 

45 _u_(<o) = a P u p (co) + a s Us(co) (2) 

The terms a p and a s are column vectors representing 
P and Sv wave particle motions at frequency co at 
each of the (2M + 1 ) measurement locations. These 
vectors include elements for M measurement loca- 

50 tions above and below the depth level correspond- 
ing to the center element in each vector. The vectors 
are defined as follows: 

a p T = (h P T d p - M ,..., h p T h p T d p M ) (3a) 

a s T = (hs T d s - M ,..., hs T ,..., hs T ds M ) (3b) 

55 where 

d p m = eticomDcosTp )/v p apc j _ ^i&imDcos^s )/v s 

are expressions in the frequency domain for the 
phase shifts of the P and S v waves at each 
measurement location relative to the center 
60 measurement location; m is the measurement 
location index -M to +M relative to the center 
location at m = 0; 

D is the spacing between measurements loca- 
tions; and 

65 h p T = (-sin *F P , cos *F P ) and h s T = (cos *F S , sin 
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M's) are unit vectors defining the x and z coordinate 
directions. 

It should be noted that the vectors a p , a s . hp and h s 
are all column vectors but are set forth herein in 
transposed form (as indicated by the T superscript) 5 
for convenient representation. Additionally express- 
ion (3a) and (3b) t and in particular the expressions 
for d p m and d s m , assume an equal spacing between 
measurement locations. This geometry is not essen- 
tial but has the advantage of simplifying the analysis. 10 

As noted above, equation (2) is not used to 
actually derive the estimated waveform components 
jj^co), because the P and Sv waves u p (co) and u s (co) 
are not directly measurable separately using geo- 
phones. Indeed, one aspect of the present invention 15 
involves a determination of these waveforms separ- 
ately from one another. However, the inventor hereof 
has found that it is possible to minimize the 
expression (1) in two steps, first for the desired 
parameters v p , 4V v s and 4V and then, if required, 20 
for the waveforms u p (co) and u s (co). 

In accordance with the invention, the first minimi- 
zation step is performed by finding values for the 
parameters v p , 4V v s and 4^ which maximize the 
value of the expression 25 

Z<o|(a P + a s ) H u(co)|2/||a p + a 3 || 2 + |(a P 
-a s ) H u(o))| 2 /||a P -a s || 2 (4) 

where () H denotes the complex conjugate transpose 
of a matrix. 

Thereafter, if desired, the P and Sv waveforms 30 
u p (o>) and u s (co) in the frequency domain may be 
separated from one another. To this end the values 
for the vectors a p and a s given by equations (3a) and 
(3b), for the values of the parameters v p , 4V v s and 
4 / s obtained from expression.^), are adopted as 35 
estimated vectors _a p and ja s . These estimated 
vectors are then incorporated in the following 
equations to obtain estimated separated waveforms 
u P jco) and u s(<o) ^ ^ ^ ^ ^ 

u p (co) = [Up + a_ s )/|U P + _a_sll 2 + U P -a_s)/lla p 40 
-a^ll 2 ] H u(co) ^ (5a) ^ ^ ^ ^ ^ 
^upfco) = [(_a P + _a s )/ll_a P .+ _a_ s ||2 -(_a p -a_ 3 )/||a p 
-_a s )ll 2 ] H u(co) (5b) 

The procedure summarized above will now be 
described in more detail with reference to Figures 4a 45 
and 4b. At step 100, in Figure 4a, a vertical seismic 
profiling operation is performed as shown in Figure 1 
to obtain waveform measurements along three 
orthogonal axes at a multiplicity of locations in the 
borehole 14. 50 

Typically the orientation of the measurement axes 
of the geophones in the tool 10 is not coincident with 
that of the coordinate system illustrated in Figure 2, 
since the orientation of the tool 10 about its 
longitudinal axis is uncontrolled. Accordingly, at step 55 
102 the waveform measurements that would have 
been obtained had the geophone axes been aligned 
with the x-z plane of Figure 2 are derived ( a 
procedure known as rotating the data). If the 
orientation of the tool 10 is monitored during data 60 
acquisition the angular offset between the two sets 
of axes is known and this offset information can be 
used in conjunction with the recorded waveforms in 
known manner to perform the rotation. If the tool 
orientation is not known it is still possible to carry out 65 



the rotation operation using information in the 
three-component measurements themselves - see 
Esmersoy. C. "Polarization analysis, rotation and 
velocity estimation in three component VSP\ in 
Vertical Seismic Profiling, Part B: Advanced Con- 
cepts, M.N. Toksoz and R.R. Stewart (Ed.). Geophy- 
sical Press, 1984. 

The (rotated) waveform data includes contribu- 
tions both from direct, down-going seismic energy 
and from reflected, up-going energy. The subse- 
quent analysis requires only the down-going compo- 
nents of the waveforms, so a velocity filter is applied 
in known manner at step 104 to separate the 
up-going and down-going components. The down- 
going components are selected for further analysis 
in accordance with this invention, while the up-going 
components may be retained for use in other 
techniques. 

Direct shear waves from the source 26 and shear 
waves resulting from conversion of compressional 
wave energy near the surface may propagate at a 
different angle 4 / s compared to shear waves result- 
ing from mode conversions at bed boundaries within 
the formation. In order to avoid obtaining a distorted 
value for the angle 4V it is therefore desired to 
separate shear wave components travelling at 
different angles and analyze them individually. 
Accordingly a technique known in the art and 
referred to as time-muting is applied at step 106 to 
the waveforms to exclude waveform components, 
such as direct shear waves, occurring outside a time 
window determined as part of the muting procedure. 

At step 108 a Fourier transform procedure is 
applied to the waveforms to convert them into the 
frequency domain for the subsequent steps of the 
analysis. 

The main procedure for evaluating v p , 4V v s and 
4 f s (steps 1 1 4 to 1 32 described below) at each depth 
level (measurement location) uses waveforms from 
a number M of depth levels above and below that 
depth level. At step 110 this number M is deter- 
mined. Typical values for M are in the range one to 
five. For a given spacing D and signal-to-noise ratio, 
a smaller value of M provides better resolution but 
the possible inaccuracy in the estimated values 
produced is larger. At step 112 the first depth level 
for which parameters are to be evaluated is selected, 
for example the depth level M + 1 measurement 
locations above the lowest such location. 

The analysis to be described below requires 
estimates of the parameters v p . 4V v s and 4 / s as a 
preliminary step to evaluation of expression (4). 
Once the derivation of these parameters for a first 
depth level has been accomplished, the derived 
values can be used as estimates for commencing 
the derivation for the next adjacent depth level. 
However, for the first depth level itself no such 
previous result exists. Accordingly at step 1 13 initial 
values appropriate to the first depth level are 
provided. Predetermined constant values of typical 
velocities and angles may be used, or specific values 
may be derived. For example, a polarization analysis 
may be made of the waveforms obtained in the 
vicinity of this first depth level, in order to derive 
initial estimates for the angles of incidence 4*p and 
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^Fs. This analysis may be performed for example as 
described in the above-mentioned Esmersoy article. 

In principle the moveout of first breaks in the 
waveforms obtained near the first depth level can 
likewise be used to obtain similar estimates for the 5 
velocities v p and v s . However, a direct computation of 
velocity, from the spacing between measurement 
locations of the tool 10 and from the moveout along 
the time axis between waveforms from adjacent 
measurement locations, will most likely provide an 10 
incorrect answer. This is because such a computa- 
tion does not allow for the possibility that the seismic 
energy is propagating in a direction other than 
parallel to the borehole 14. Thus the estimates *F P 
and 4*5 of the angles of the actual directions of 15 
propagation must be incorporated in any initial 
estimation of the wave velocities from the first break 
moveouts. 

The procedure then commences evaluating the 
required parameters for the first depth level, enter- 20 
ing a loop at step 114. At this step the estimated 
values for v p , Tp, v s and 4^ used to evaluate the 
factors d p m and d s m appearing in expressions (3a) 
and (3b) as well as the vectors h p and h s , and these 
factors and vectors are then used in evaluating 25 
expressions (3a) and (3b) to obtain the vectors a p 
and a s . As indicated above, upon the first pass 
through the loop starting at step 114 the estimated 
values used are those obtained at step 113. 

At step 1 1 6 the waveform measurements for the M 30 
measurement locations above and below the pres- 
ently selected depth level are selected to form for 
each frequency co the complex column vector u(co). 
This vector and the vectors a p and a s determined at 
step 114 are then applied to the evaluation of 35 
expression (4), at step 118. In this connection it is 
useful to note that the denominators of the two 
fractions in expression (4) can be directly evaluated 
from the estimated angles and as follows : 

|[a P + a s |p = 2(2M + 1) + 2sin(¥ s - x F p )sin[(2M + 40 
1 )a]/sinoc (6a) 

;!a P -as|| 2 - 2(2M + 1) -2sin(*F s - 4 / p )sin[(2M + 
1)a]sina (6b) 

where a 1 / 2 coD(cosHVv s - cos*F P /v p ). 
Evaluation of the numerators of the fractions in 45 
expression (4) for a single frequency co involves 
adding or subtracting the column vectors a p and a s 
element by element, changing the sign of the 
imaginary part of each element in the resulting sum 
or difference vector to produce its complex con- 50 
jugate and transposing the elements to produce a 
row vector, multiplying this row vector by the column 
vector u(co) for the frequency co and obtaining the 
square of the magnitude of the resulting scalar 
quantity. Dividing the squared magnitudes by the 55 
appropriate denominator value from (6a) or (6b) 
gives the values of the two fractions which are added 
and then summed with the corresponding values for 
other frequencies co to produce a value for express- * 
ion (4). 60 

A test is performed at step 120 to determine 
whether the maximum value of expression (4) has 
been found. This may be done, for example, by 
comparing the value obtained at step 118 with the 
value obtained in previous passes through the loop 65 



(except of course for the very first pass). If the value 
has changed by an amount less than some predeter- 
mined threshold (for example by less than 10 -5 of 
the previous value), then the value obtained at step 
118 is accepted as the maximum value, and the 
estimates of the parameters v p , *F Pl v s and *F S used to 
obtain that value are accepted as the actual values of 
those parameters. 

If the test at step 120 does not find a sufficiently 
small change in the value of expression (4), the 
procedure returns to step 114 via step 122 where 
new estimates for the the parameters v p , *F Pl v s and 
*F S are derived. These estimates may be derived 
using known techniques such as a quasi-Newton 
method or the Fletcher-Poweli method (Fletcher, R., 
"Practical methods of optimization - Volume 1, Un- 
constrained Optimization", John Wiley and Sons 
1980). Such techniques typically involve evaluating 
the (partial) derivatives of expression (4) with 
respect to each of the parameters v p , 4^, v s and ^Fs 
and testing them to find a combination of changes of 
appropriate magnitude in the values of these 
parameters which will alter the value of expression 
(4) towards its maximum. 

When the test at step 120 indicates that the 
maximum value of expression (4) has been obtained 
to the desired degree of accuracy, the procedure 
continues to step 124 where the corresponding 
estimates of v p , 4^, v s and are recorded as the 
actual values of these parameters. Thereafter the 
values of the vectors a p and a s evaluated at step 1 14 
for these parameter vajues are^used at step 126 as 
the estimated vectors _a p and ja_ s to derive separate 
compressional and shear waveforms u p (co) in 
accordance with equations (5a) and (5b). These 
separate waveforms are likewise recorded for 
further use. 

At step 128 a test is performed to determine 
whether the waveforms for all depth levels have 
been analyzed. If not the procedure advances to 
step 1 30 where the waveform data for the next depth 
level is brought in to the group of 2M + 1 waveform 
sets, replacing the data for the depth level at the 
other end of the 2M + 1 measurement locations. 
The values for the parameters v Pl W Pt v s and 4*5 found 
at step 124 are adopted as first estimates of their 
values at the next depth level, at step 132, and the 
procedure then returns to step 114 to evaluate the 
parameters. 

When the test at step 128 determines that the data 
for all depth levels has been analyzed the procedure 
terminates at 134. 

There has been described and illustrated herein a 
method in accordance with the present invention for 
evaluating parameters related to the elastic proper- 
ties of subsurface earth formations. While particular 
embodiments of the invention have been described, 
it is not intended that the invention be limited 
thereby. Thus the description above envisages the 
measurement of seismic energy along three ortho- 
gonal axes with the tool 10. However, it is possible to 
practice the invention in the analysis of measure- 
ments obtained along only two orthogonal axes, if 
these axes are aligned with the x and z directions as 
shown in Figure 2. Therefore it will be apparent to 
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those skilled in the art that various changes and 
modifications may be made to the invention as 
described without departing from the spirit and 
scope of the appended claims. 



Claims 



10 

1. A method of evaluating parameters related 
to the elastic properties of a subsurface earth 
formation traversed by a borehole, comprising 
the steps of: 

obtaining plural-axis measurements from a 15 
seismic profile operation in said borehole; 

describing seismic propagation through 
said formation in terms of a model incorporating 
a plane compressional wave and a plane shear 
wave propagating with unspecified velocity and 20 
direction; 

deriving values for the velocity and direction 
of propagation of said plane waves that provide 
a best fit of said model to said measurements; 
and 25 

producing a tangible record of said derived 
values. 

2. The method of claim 1 , wherein said values 
are derived to provide a best fit of said model to 

the measurements obtained at a small number 30 
of adjacent measurement locations. 

3. The method of claim 1 , wherein said values 
of compressional wave velocity v p and angle of 
incidence *F P at the borehole and of shear wave 
velocity v s and angle of incidence *F S at the 35 
borehole are derived in accordance with the 
expression: 

maXvp.Tp.vs.H'sSojKap + a s ) H u(co)| 2 /lla p + 
a s ||2 + |(ap-a s ) H u(o>)| 2 /lla p -asl|2 

where 40 

co is angular frequency; 

a p and as are column vectors representing 
compressional and shear wave particle motions 
at frequency co at a multiplicity of locations in 
said borehole; 45 

u(co) represents components in a plane of 
waveform measurements obtained at said loca- 
tions; 

(i || represents the Euclidean norm of a 
complex vector; and 50 

() H denotes the complex conjugate trans- 
pose of a matrix. 

4. The method of claim 1, wherein said 
measurements are obtained along three ortho- 
gonal axes. 55 

5. The method of claim 1 , including the step of 
separating components of said measurements 
corresponding to shear waves propagating in 
different directions, before said values are 
derived. 60 

6. The method of claim 3, wherein said 
measurements are obtained at locations equally 
spaced in said borehole and said column 
vectors a P and a s are evaluated in accordance 

with the expressions 55 



a p T = (h D T dp- M h p T h p T dp M ) 

a s T = (hs T d s ' M h s T h s T d s M ) 

where 

d p m = e l«»mDcosM'p j v p and (J^ = e |,,: ' mD 

cosf s i-vs are expressions in the frequency 
domain for the phase shifts of compressional 
and shear waves at each measurement location 
relative to a center measurement location ; 

m is a measurement location index -M to 
-f M relative to the center location at m = 0; 

D is the spacing between measurement 
locations and 

h p T « (-sin T p , cos H'p) and h s T = (cos Ts, 
sin ^s) are unit vectors defining two orthogonal 
coordinate directions. 

7. The method of claim 1 , including the step of 
deriving separate compressional and shear 
waveforms corresponding to the values derived 
for the velocity and direction of propagation of 
said plane waves. 

8. The method of claim 7, wherein said 
separate compressional and shear waveforms 
u p (co) and u s (co) are derived in accordance with 
the expressions 

^ up^co) = ((_a P + _a s )/IU P + _a*ll 2 + (a_ p 
-±s) y /lla_p-_a s ||2j^u(0)) ^ ^ 

u P £io) = [(a_ p + j_ s )/ll_a_p + _a s ll 2 -(j_ P 
-±s)/|l^-JLsll 2 ] H u(co) 

where a p and a s are values of the column 
vectors a p and a s corresponding to the values 
derived for the velocity and direction of propa- 
gation of said plane waves. 

9. The method of claim 1 , including the step of 
deriving estimates of the direction of propaga- 
tion of said plane waves by a polarization 
analysis of said measurements. 

10. The method of claim 9, including the step of 
deriving initial estimates of the velocity of said 
plane waves by analysis of moveout in said 
measurements and of said initial estimates of 
the direction of propagation of said plane 
waves. 
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Obtain 3-axis waveforms 
from vertical seismic 
profiling operation 
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Rotate waveform data to 
derive components in x-z 
plane 



Ve 1 oc i ty filter 
waveforms to separate 
up-going (reflected) and 
down-going components; 

select down-going 
components for analysis 



Fig. 4a 
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Apply time-muting to 
eliminate direct shear 
wave components 
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Fourier transform al 1 
waveform records 
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Choose span -M to of 

waveforms to be 
incorporated in 'analysis 
for each depth level 
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Select starting oeDZn 
level: M+l measurement 
locations from lowest 
level 
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Provide initial values 
for Vp, iif D , v s and u/ s 
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Use estimated values of 

v p , ip p . v s and ip s to 
evaluate a p and a s from 
expressions (3a) & (3b) 

— + 
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Select waveform 
measurements from -M to 
+M relative to present 
depth level to form u(co) 

T 
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Evaluate expression (4) 
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Derive new estimated 
values v p , \|i p , v s and u/ s 
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Record values or v p . u/p . 
v s and u/s 
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Derive separate u 0 (co) 
and u s (co) waveforms from 
expressions (5a) & (5b); 
record waveforms 




Shift in waveform data 
for next depth level 
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Fig.4b 
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Adoot values of Vr 



v s and U/ s from steo i24 
as new estimated values 
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@ Three-axis vertical seismic profile measurements 
are analyzed to derive values for the velocity v p and 
angle of incidence * p of a plane compressional wave 
at the borehole and for the shear wave velocity v s 
and angle of incidence * s of a plane shear wave at 
the borehole. In the described embodiment this is 
accomplished by finding values for these parameters 
which maximize the expression 



EJ( a P 
(«)| 2 Ha p -a JF 
where 



+ a s ) H u M| 2 /lla p 
(4) 



+ a s ll 2 + |( a p -a s ) H u 



w is angular frequency; 

a p and a s are column vectors representing 
compressional and shear wave particle motions at 
frequency u at a multiplicity of locations in the 
borehole; 

u (cj) represents components in a plane of 
waveform measurements obtained at these locations; 

II II represents the Euclidean norm of a complex 
vector; and 

() H denotes the complex conjugate transpose of a 
matrix. 
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